Stationary Scalar Clouds Around Rotating Black Holes 
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Motivated by novel results in the theory of wave dynamics in black-hole spacetimes, we analyze 
the dynamics of a massive scalar field surrounding a rapidly rotating Kerr black hole. In particular, 
we report on the existence of stationary [infinitely long-lived) regular field configurations in the 
background of maximally rotating black holes. The effective height of these scalar "clouds" above 
the central black hole is determined analytically. Our results support the possible existence of 
stationary scalar field dark matter distributions surrounding rapidly rotating black holes. 



I. Introduction. The 'no-hair' conjecture [3, 0] 
has played a key role in the development of black-hole 
physics UQ. This conjecture asserts that black holes 
are fundamental objects — they should be described by 
only three externally observable (conserved) parameters: 
mass, charge, and angular momentum. 

The physical idea behind the no-hair conjecture is 
based on an intuitive (and sometimes oversimplified) 
picture according to which all matter fields which are 
present in the exterior of a black hole would eventu- 
ally be swallowed by the black hole or be radiated away 
to infinity (with the exception of fields which are asso- 
ciated with conserved charges) @, 0- In accord with 
this logic, various no-hair theorems indeed exclude static 
(time- independent) scalar fields massive vector fields 
0, and spinor fields 10] from the exterior of black holes. 

It should be stressed, however, that the various no-hair 
theorems [ll-flo| do not rule out the existence of time- 
dependent field configurations in the black-hole exterior. 
Indeed, in a very nice work Barranco et. al. llj have 
recently found time- decaying regular scalar field configu- 
rations surrounding a Schwarzschild black hole that can 
survive for relatively long times (as compared to the dy- 
namical timescale set by the mass of the black hole) . 

Ultra-light scalar fields have been invoked in recent 
years as possible candidates to play the role of the dark 
matter component of our universe (see llMl4| and refer- 
ences therein) . Given the fact that most galaxies seem to 
contain a super-massive black hole at their centers [lij ]. 
it was pointed out in ll| that in order to be a viable 
candidate for the dark matter halo a scalar field config- 
uration must be able to survive for (at least) cosmologi- 
cal timescales. In this context, the important conclusion 



presented in 11| is very encouraging: time-decaying con- 
figurations made of ultra-light scalar fields surrounding a 
supermassive Schwarzschild black hole can indeed survive 
for such extremely long periods. 

The main aim of the present work is to extend the im- 



portant results of ll|_m three new directions: 
(1) The analysis of llj was restricted to a spherically 
symmetric Schwarzschild black-hole background. How- 
ever, it is well-known that realistic black holes generally 



rotate about their axis and are therefore not spherical. 
Thus, an astrophysically realistic model of wave dynam- 
ics in black-hole spacetimes must involve a now-s phe rical 
background geometry with angular momentum [16j |. In 
the present work we consider such non-spherical back- 
grounds and analyze the dynamics of massive scalar fields 
in realistic rotating (Kerr) black-hole spacetimes. 



(2) The regular scalar configurations considered in [11 1 
have long yet finite lifetimes. In the present work we 
shall prove the existence of stationary (i.e. infinitely long- 
lived) regular field configurations surrounding rotating 
black holes. 



(3) The focus in [llj was mainly on black hole-scalar field 
configurations characterized by the dimensionless prod- 
uct Mfi <C 1, where M is the mass of the central black 
hole and /j, is the mass of the surrounding scalar field. 
In the present work we shall analyze the complementary 
regime of M/i > i. 

II. Description of the system. The physical sys- 
tem we consider consists of a massive scalar field coupled 
to a rotating Kerr black hole of mass M and angular mo- 
mentum per unit mass a. In order to facilitate a fully 
analytical study, we shall assume that the black hole is 
maximally rotating with a = M. Similar results can be 
obtained (with the cost of a more involved analysis) for 
non-extremal rotating black holes. In the present study 
we shall restrict ourselves to the case of a test scalar field 
in the background of the rotating Kerr black-hol e sp ace- 
time. We shall therefore use the terminology of [llj and 
talk about scalar "clouds" surrounding the black hole 
rather than a genuine scalar "hair" 
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The dynamics of a massive scalar field ^> in the Kerr 
spacetime is governed by the Klein-Gordon equation [lil ] 



(V a V a - M 2 )* = 



(1) 



(Here /i stands for A4G/hc, where M. is the mass of the 
scalar field. We shall use units in which G = c = % = 1. 
In these units \i has the dimensions of 1/length.) One 
may decompose the field as 



* lm [t,r,9,4>) = e" n *Si m (0; Mcj)R lm (r ; M, u)e 



(2) 
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where (t, r, 9, <p) are the Boyer-Lindquist coordinates, ui is 
the (conserved) frequency of the mode, I is the spheroidal 
harmonic index, and m is the azimuthal harmonic index 
with — I < m < I. (We shall henceforth omit the indices I 
and m for brevity.) With the decomposition @, R and 5 
obey radial and angular equations both of confluent Heun 
type coupled by a separation constant K(Mto) [lil-l^. 

The angular functions S(9; acu) are the spheroidal har- 
monics which are solutions of the angular equation 
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The angular functions are required to be regular at the 
poles 9 = and 9 = 7T. These boundary conditions 
pick out a discrete set of eigenvalues {Ki m } labeled by 
the integers I and m. For M 2 (/x 2 — a; 2 ) 



< 



m one can 



treat M 2 (cj 2 — [i 2 )cos 2 9 in Eq. (J3|) as a perturbation 
term on the generalized Le gen dre equation and obtain 
the perturbation expansion [21 [ 



oo 



Ki m +M (fj, —ui ) = 1(1 + 1)+) c k M 2k (^-u z Y (4) 



fe=i 



for the separation constants K\ m . The expansion coeffi- 
cients {ck(l,m)} are given in Ref. (2lj . 



The radial Teukolsky equation is given by [lg, |22j, [23 1 



c?r ( dr ) 



# 2 +A[2mAfw-A-/i 2 (r 2 +M 2 )] 



R 



= 0, 
(5) 

where A = (r — M) 2 and = (r 2 + Af 2 )w - mM. The 
degenerate zero of A, rn = M, is the location of the 
black-hole horizon. 

We are interested in solutions of the radial equation 
([5]) with the physical boundary conditions of purely in- 
going waves at the black-hole horizon (as measured by 
a comoving observer) and a decaying (bounded) solution 
at spatial infinity [24j-[29(. That is, 



as r — > oo (y — > oo) ; 



^ e -*("-"*n)» asr ^ rff (y^-oo), 

where the "tortoise" radial coordinate y is defined by 
dy = [(r 2 + M 2 )/A]dr. Here Q = 1/2M is the angular 
velocity of the black-hole horizon. 

Note that a bound state (a state decaying exponen- 
tially at spatial infinity) is characterized by to 2 < [i 2 . For 
a given mass parameter /i, the boundary conditions (0 
single out a discrete set of resonances {w„(/i)} which cor- 
respond to the bound states of the massive field [US!. 
Stationary resonances, which are the solutions we are in- 
terested in in this paper, are characterized by Quj = 0. 
(We note that, in addition to the bound states of the 



massive field, the field also has an infinite set of discrete 
quasinormal resonances 33~33j which are characterized 
by outgoing waves at spatial infinity.) 

III. The stationary scalar resonances. As we 
shall now show, the field @ with u = mil describes a 
stationary regular solution of the wave equation (fT]). It 
is convenient to define a new dimensionless variable 



r - r H 
x = , 

in terms of which the radial equation ([5]) becomes 

,d 2 R „ dR Tr „ „ 
x 2 — + 2x— + VR = , 

dx z dx 



(7) 



(8) 



where V = (m 2 j A-M 2 [?)x 2 + (m 2 -2M 2 p?)x+(-K lm + 
2m 2 -2M 2 [i 2 ). Defining 



f = xR and z = 2y//j? -m 2 /AMx , (9) 
one obtains the radial equation 



dz 2 



1 K 



f = 



(10) 



with 



TO 2 - 2M 2 il 2 , „o r , 1 „ 2 „ , ,2 2 

K = = — and 8 = K hn -\ 2m 2 +2M 2 n 2 . 

^AM 2 ^ 2 - m 2 4 

(11) 

We shall assume without loss of generality that > 0. 
Equation (TlOl) is the familiar Whittaker equation; its so- 
lutions can be expressed in terms of the confluent hyper- 
geometric functions M(a,b,z) [21, 28, 34 1 



R = Az-^ +f3 e-^ z M(- + /3-K ) l + 2@,z) + B(l3 -> -8) , 

(12) 

where A and B are constants. The notation (8 — > —0) 
means "replace 8 by —8 in the preceding term." 

The near-horizon (z -+ 0) limit of Eq. (|12l) yields 



21|,[34| 



R 



Az-i +f> + Bz-v-P 



(13) 



Regularity of the solution at the horizon (R is bounded 
for z —¥ 0) requires 



B = and ®.B > 



1 



(14) 



Approximating Eq. (TT21 for z — > oo one gets [2l|, \34 
T(1 + 2B) 



R^ A 

r(i + 2/3) _ 
'ra + 8 + K ) z 



ra + 8-n) 



(15) 



A bound state is characterized by a decaying field at spa- 
tial infinity. The coefficient 1/T(^ + 6 — n) of the growing 
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exponent e' 2 in Eq. (1151) should therefore vanish. Using 
the well-known pole structure of the Gamma functions 
21] . we find the resonance condition for the stationary 
bound-states of the field: 



(16) 



where the resonance parameter n is a non-negative inte- 
ger (n = 0, 1, 2, ...). Substituting Eqs. ^ and ^ into 
Eq. (TT2")) . one obtains the compact form 



R = Az-^e-^L^\z) 



(17) 



for the radial solutions, where Ln^\z) are the general- 
ized Laguerre Polynomials (see Eq. 13.6.9 of [21(). 

We note that the r.h.s of the resonance condition (fl6|) 
is positive definite. Taking cognizance of Eq. (|lip with 
k > 0, one concludes that the bound-state resonances 
must lie within the band 



to , , to 
-2 <M »<T2 



(18) 



Note that the inequality (fl8|) excludes the existence of 
regular static solutions with m = (lj = 0). Of course, 
this finding is in agreement with the well-known no-hair 
theorems which indeed exclude static hairy configu- 
rations. However, as we shall now show, time-dependent 
(stationary) resonances do exist! 

In order to solve the resonance condition, Eq. (p~6|) . it 
is instructive the introduce the dimensionless variable 



e = v/A/V - (m/2) 2 , (19) 
in terms of which k and /3 can be expressed as 
(to/2) 2 - 



and /3 Z 



<'+5> 2 - 



' 2,2 
-TO + C - 



-E 

fc=i 



Cfe£ 
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(20) 

Substituting (f2TI|) into (flB]) , one finds that the resonance 
condition can be expressed as a polynomial equation for 
the dimensionless variable e 35]: 



?7i 4 -4TO 2 (2n+l)e + 16 



TO 2 -(; + i) 2 + (n + i) 2 



+16(2n + l)e 3 - 16cie 4 - 16 ^ 



2fc+2 _ 



0(21) 



k=2 



To demonstrate the existence of a discrete and infinite 
family of stationary resonances, we present in Table U 
the dimensionless quantity Mfj, for the fundamental fam- 
ily (I = m = 1) of resonances. We displa y r esults for 
various values of the resonance parameter n [36| . A qual- 
itatively similar behavior is observed for other families of 
resonances (characterized by the two parameters {I, to}). 

The infinite spectrum {Mfi(n)} of field-masses satisfy- 
ing the resonance condition (1161) is a decreasing function 



(for n > 1) of the resonance number n. For n 3> I one 
finds from (|2T1) 



to 

y 



for the resonances 137 



16n 2 8n 3 



+ 0(n~ 4 ) 



(22) 



n 


resonance 


^cloud 





0.526 


8.557 


1 


0.510 


23.485 


2 


0.505 


46.398 


3 


0.503 


77.313 



TABLE I: Stationary scalar resonances of a maximally rotat- 
ing Kerr black hole. We display the dimensionless product 
M/j, for the fundamental family of resonances (I = m = 1) 
and for various values of the resonance parameter n. Also 
shown are the effective heights (in units of the black-hole ra- 
dius), iEcioud, of the scalar clouds above the central black hole. 
In general, the size of a scalar cloud increases with increasing 
resonance number n. Note that all resonances conform to the 
lower bound (l23l). 



IV. Effective heights of the scalar clouds. We 

shall next evaluate the effective heights of the stationary 
scalar "clouds" which surround the central rotating black 
hole. These clouds correspond to the family of wave- 
functions ([T7"l) satisfying the resonance condition (fl6|) . 

It is worth mentioning that a nice 'no short hair theo- 
rem' was proved in 0, |38| for static and spherically sym- 
metric hairy black-hole configurations. According to this 
theorem, the "hairosphere" [the region where the non- 
trivial (non-asymptotic) behavior of the black-hole hair is 
present] must extend beyond | the horizon radius. This 
translates into the lower bound 



1 

> - 

~ 2 



(23) 



where x is the dimensionless height defined in (|7J). A 
rough estimate of the size of our stationary scalar config- 
urations can be obtained by defining their effective radii 
as the radii at which the quantity 47rr ,2 |4'| 2 attains its 
global maximum. Taking cognizance of Eq. (|17j) with 
Lq\z) = 1, one finds 



c (0) 

cloud 



2,3+1 
2e 



(24) 



for the ground-state resonances with n — 39]. For 
higher n values one generally obtains larger scalar clouds 
with the approximated relation a^ ud ~ 2p+i+2n ^ gee 

Eq. erg. 

The effective heights of the fundamental scalar clouds 
above the rotating central black hole are given in Ta- 
ble HI from which we learn that {x c \ ou d(n)} are always 
larger than i. For n » / one finds from (f2"2")l a^ioud — 
4n 2 /m 2 ^> 1. It is remarkable that our stationary scalar 
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clouds respect the lower bound (|23j) despite the fact that 
they do not satisfy the conditions of the original theorem 
@: they are neither static nor spherically symmetric. 

V. Summary and discussion. In this Letter we 
have analyzed the dynamics of a massive scalar field in an 
astrophysically realistic (rotating) Kerr black- hole space- 
time. In particular, we have proved the existence of a 
discrete and infinite family of resonances describing sta- 
tionary (non-decaying) scalar configurations surrounding 
maximally rotating black holes. The effective heights of 
these scalar "clouds" above the central black hole were 
determined analytically - it was shown that these non- 
static and non-spherically symmetric configurations con- 
form to the lower bound (|23|) originally derived in [6| for 
static and spherically symmetric hairy black-hole con- 
figurations. Thus, our analysis provides the first direct 
evidence for the general validity of the bound (|2T>1) . 

Our results support the possible existence of station- 
ary scalar field dark matter distributions surrounding 
astrophysically realistic (rota ting ) black holes. More- 
over, while former studies [ill 0] focused on the regime 
A//i <C 1, our analysis opens the possibility for the ex- 
istence of such black hole-scalar field configurations in 
the new regime M/i > 1/2. Our analytical findings for 
maximally rota ting black holes are in accord with the nu- 
merical work of [27] who discussed slowly-decaying bound 
states of non-extremal black holes. 

Probably the most interesting question is the following 
one: what is the final end-state of gravitational collapse? 
The no-hair conjecture [l[ asserts that the final outcome 
of rotating gravitational collapse is a bald Kerr black 
hole. However, our results suggest the following alter- 
native scenario: during the early stages of the evolution 
some of the fields would indeed be radiated away to in- 
finity while some would be swallowed by the newly-born 
black hole (see ll| for the spherically symmetric case). 
However, if the initial matter distribution contains mas- 
sive scalar fields, then those fields which satisfy the res- 
onance condition (j!6[) may actually survive as infinitely 
long-lived (stationary) resonances outside the black hole. 
Thus, the final configuration is expected to be a rotat- 
ing black hole surrounded by stationary scalar clouds. It 
would be interesting to verify t his p rediction using nu- 
merical simulations as done in 11] for the spherically 
symmetric (non-rotating) case. 
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